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Abstract. The two-nucleon contributions to the electric dipole moment (EDM) of the deuteron, induced 
by the QCD #-term, are calculated in the framework of effective field theory up-to-and-including next- 
to-next-to-leading order. In particular we find for the difference of the deuteron EDM and the sum of 
proton and neutron EDM induced by the QCD #-term a value of ( — 6 ± 4) 9 x 10~ 4 efm. The by far 
dominant uncertainty comes from the CP- and isospin-violating nNN coupling constant. Matrix elements 
are evaluated with an accuracy of better than 10%. 
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Electric and magnetic moments - 24.80. +y Nuclear tests of fundamental interactions and symmetries - 
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1 Introduction 

Under the assumption that the CPT theorem is valid, per- 
manent electric dipole moments (EDMs) of elementary 
particles and nuclei, which arise under parity P and time- 
reflection T breaking, belong to the most promising signals 
of CP-violating physics beyond the Cabibbo-Kobayashi- 
Maskawa (CKM) phase of the Standard Model (SM) [TJH 
[3] . Possible mechanisms [HG2 are the dimension- four 8 vac- 
uum angle term of Quantum Chromodynamics (QCD) [5] 
and the effective dimension-six quark, quark-color, and 
gluon-color terms [TUHJH] (including certain combinations 
of four-quark terms [10] ) resulting from extensions of the 
SM such as supersymmetry |llj , many-Higgs scenarios |12j 
etc. In refs. [T31H¥| it was recently pointed out that the 
same mechanism that drives the potential CP violation 
beyond the SM in D -> K+K-fir+n- should, if 

present, also lead to an enhanced nucleon EDM. How- 
ever, a single successful measurement of an EDM signal of 
the neutron, say, would not suffice to isolate the specific 
CP-violating mechanism. Therefore, more than one EDM 
measurement involving other hadrons and (light) nuclei, 
e.g. the proton, deuteron, helium-3, are necessary in order 
to uncover the source(s) of the CP breaking. 

In recent years various theoretical studies focussed on 
the calculation of EDMs for light nuclei [17IU8UTgp O PH 
22.23,24 , largely triggered by on-going plans for dedi- 
cated experiments to measure EDMs of light ions using 
storage rings 25,26,27,28,29 . These calculations revealed 



that different CP-violating mechanisms contribute to dif- 
ferent probes with different strength. Therefore, non-zero 
measurements as well as controlled calculations of nucleon 
and nuclear EDMs are necessary to reveal additional in- 
formation on the physics beyond the SM that drive non- 
vanishing EDMs. 

In this work we calculate the two-nucleon contribu- 
tion to the deuteron EDM that would be produced from a 
non-vanishing QCD 6>-term up to next-to-next-to-leading 
order. Thus, once the EDMs of the proton, neutron and 
deuteron were measured, the results of our calculation 
would allow one to extract the value of 6 directly from 
data, assuming that no other CP-violating mechanisms 
contribute significantly. Since lattice QCD will eventually 
be able to calculate the neutron and proton EDMs with 
6 as the only input, a combination of the calculation pre- 
sented here with lattice QCD and experimental numbers 
will enable one to decide, if the 9-term is the culprit of gen- 
erating the EDMs. Note that direct lattice calculations for 
nuclear EDMs would be much more challenging. 

The terms of the CP-violating interaction Lagrangian 
relevant for this work are given by — see ref. |22j and 
references therein — 

= N (b + &1T3) S^Nv'F^ + 9o Ntv ■ tN + 9l N^N 
+C°NND fl (NS> 1 N) + C^NtND^NtS^N) 
+CfNT 3 ND fl (NS f *N) + C$NND fl (NT 3 S' i N) 
+ .... (1) 
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Fig. 1. Total CP-violating transition current. The P- and T-violation stem from either CP-violating two-nucleon potentials or 
two-nucleon-irreducible CP-violating transition currents. 



Here v 11 = (1,0), = (0, |cr) and r are the nucleon 
velocity, spin and isospin, respectively, while is the co- 
variant derivative. For #-term-induced CP violation naive 
dimensional analysis gives that g\jg\ ~ eM 2 /m 2 N [55] . 
However, as already pointed out in ref. [18) and refined 
further below, gf is significantly enhanced compared to 
this estimate — in fact, the contribution from gf domi- 
nates the deuteron EDM. 

The single-nucleon EDM from the (9-term starts to con- 
tribute at the one loop level [3J[3]. At the same order 
there are two counter terms, proportional to bo and b\, 
to absorb the divergence [50115X1152"] . Therefore, although 
the value of the CP-violating coupling constants go and 
gi can be related to the strength of the QCD 6>-term, 9, 
within the effective field theory the same is not possible for 
the EDM of a single neutron or proton. This is different 
in case of the nuclear EDMs: for the few-nucleon contri- 
butions, counter terms appear only at sublcading orders 
and therefore controlled calculations become feasible, al- 
though, in case of the deuteron EDM, with a sizable uncer- 
tainty. Such sublcading terms can be found in the second 
and third lines of eq. (UJ), where the two terms in third 
line arc additionally suppressed by isospin breaking. The 
dots in cq. (jTJ) denote further CP-violating terms that do 
not contribute to the deuteron EDM at orders considered 
in this work. These terms include CP-violating NNttj-, 
NNittt-, NNttttj-, 4iV7-terms, CP-violating photon-two- 
pion terms, and CP-violating pure pion-terms (see ref. |33j 
for the latter class). 

There are two types of contributions that are relevant 
for the present study, namely CP-violating TV TV interac- 
tions and CP-violating irreducible TV TV — > TVTV7 transition 
currents — c.f. fig. [TJ As will be outlined below, for the 
deuteron EDM the latter kind of contributions contains 
at its leading non-vanishing order loop diagrams that are 
calculated in this work for the first time. 

The paper is structured as follows: in sect. [5] the prcf- 
actors of the CP-violating nNN couplings g$ and g\ are 
derived from the QCD 0-term. After this, a brief discus- 
sion of the power counting is presented in sect. [3] Sec- 
tion^?] contains the derivation of the two-nucleon contribu- 
tions to the deuteron EDM induced by the 0-term, where 
the TV TV potential and transition current contributions are 
discussed in subsections 14.11 and 14.21 respectively. Sec- 
tion [5j contains a brief review about the changes if the 
the 0-tcrm contributions are replaced by dimension-six 



quark-color operators. Finally, in sect.|5]a short summary 
of the presented results and an outlook are given. Ap- 
pendix [A] presents two further alternatives to derive the 
CP-violating coupling constant gf , an update of the origi- 
nal derivation by Lebedev et al. |18| and a derivation in the 
framework of SU(3) chiral perturbation theory (ChPT), 
similarly to the one of g e by [5MJ1I5T1 I52] . 



2 CP-violating tzNN couplings from the 6 
term 

On the quark level the effect of the 0-term can be writ- 
ten as mj)qi^q [5], with the reduced quark mass m» = 
m u m d /(m u +m d ) = (m u +m d )(l-e 2 )/4, where e = (m u - 
m,d)/{m u + md) = —0.35 ± 0.10 [33]. It thus behaves un- 
der chiral rotations identically to the quark mass term and 
can be included in the chiral Lagrangian via 



X± 



u*x u i ux^u with x ~ 273 (s + ip) , (2) 



where s may for our purposes be identified with the quark 
mass matrix, which reads 



M 



m d 



m d 



-T3 



(3) 



while p — m^QXi. The pion fields are contained in the 
usual SU(2) matrix u = J7 1 / 2 , see e.g. [35]. In our work 
we extend its standard representation by an additional 
pscudoscalar-isoscalar field, fj, to an U(2) matrix, such 
that it reads in the exponential representation 



exp[i'(7r-T + 77l 2 )/(2i^)] . 



(4) 



The fj field is the U(2) remnant of the U(3) rj and 7/ or of 
their U(2) component r\ n ~ (uu + dd)/\^2 (see e.g. [36,37, 
[35]). For the mass of the latter, driven by the anomaljQ, 
one finds 



M 2 



-M 2 + -M 2 , 
3 " 3 " 



(711.2 MeV) 5 



(5) 



Alternatively to the methods used in ref. [33] , the lead- 
ing CP-violating coupling constants can be derived from 



1 Thus, in a complete treatment an explicit mass term for r\ n 
should be added to the Lagrangian. 
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the ci- and the cs-tcrm of [35| . namely 



Cl N\ X +)N + c 5 N^ 
— > a AB 

+ c 5 2B TVt 



x+ - y(x^ 



N 



, 2m*0 _ 
(m u +m d ) H = — V 



N 



2m*6 

{m u -m d )T 3 H — (7T • t) 



Here (•) denotes the trace in flavor space. The dots indicate 
that terms not relevant for this study were omitted. 

We start with a discussion of the term proportional 
to C5. The first term of the third line of eq. (JSJ leads to 
the quark-mass-induccd part of the proton-neutron mass 
difference, 5m^p S . It can be quantified from three differ- 
ent sources: (i) the use of dispersion theory to quantify 
the electromagnetic part of the proton-neutron mass dif- 
ference [3HDUH1II2], (ii) lattice QCD or (iii) from 
charge-symmetry-brcaking (CSB) studies of pn —> dir° [44] , 
All analyses lead to consistent results, with the first one 
being the most accurate. Thus we will use [55] 



4B(m u - m d )c 5 = 5m s ^ g 
From this we get 



(2.6±0.5)MeV . (7) 



9o 



4FV 



where we used F v = 92.2 MeV [34]. The superscript 9 in- 
dicates that we here only include the strength that comes 
from the 0-term. The expression given above agrees with 
the prediction of rcf. [45] when eq. (14) of ref. [45] is in- 
serted into the corresponding eq. (80- It turns out that 
the value of g® is more than a factor of 10 smaller than 
the estimate from the naive dimensional analysis given by 
0M% I (ttimF^) in terms of the pion mass M T , the nucleon 
mass mjv and the pion axial decay constant F^. 

The first term in the second line of eq. © leads to the 
quark-mass-induced isoscalar contribution to the nucleon 
mass — thus c\ can be related to the irN sigma-term. For 
this low-energy constant we use the value given in ref. |46j , 



ci = (-1.0 ±0.3) GeV" 



(9) 



which is a compilation of various extractions of c\ [17,48, 
49,50]. Inspired by ref. [TB] we will in the following relate 
the amplitude of eq. ([6]) to an effective CP-violating irNN 
vertex, see fig. [U In order to do so, we still need to relate 
the fj field to the pion field. The formalism outlined above 
provides for the 7r°-ry mixing amplitude 



B(m d - mj 
Ml 



Ml 



M? 



M? 



Ml 



(10) 



2 Note that the result of ref. [35] has the opposite sign to ours 
(which is compensated by the opposite sign of e). Furthermore, 
F w is defined twice as large there. 



I ^0 
I 



! ^ 



TV Fig. 2. CP-violating tt°NN vertex g 1 (black square) in- 
duced by -K°-f)- mixing and the CP-violating fjNN vertex (open 
(6) square). 



where we used M% = 2Brh+0(M 2 ) with™ = ^(m u +m d ). 
Then we get from combining these ingredients 



9l 



Mi 



Ml 



Ml 



6 = (0.0035±0.0018)<9. (11) 



In particular, we find for the ratio 



'A 

9 6 o 



8 Cl e 2 M* 



6 (M?-Ml)6mZ 6 



-0.20 ±0.13. 



(12) 



(-0.018 ± 0.007) 9 , (8) gO/gf. 



The relatively large uncertainty of the given quantities 
arises, since we implicitly employed SU(3) ChPT in order 
to fix g\. In addition, there is some freedom in choosing 
the mass for the fields f) and r\ n : the uncertainty given 
includes Mfj = M Vn as well as Mfj = M 7 y. The value of 



<?i/<?o i s numerically about a factor of 25 larger than the 
SU(2) estimate of order eM^/m^, which is a consequence 
of the first relation of eq. p^]) if the scaling 8m s ^ ~ eM 2 

is assumed. The main origin of this difference is that g$ 
is unusually small — instead of two powers in the count- 
ing the relative suppression numerically is of the order 
of one power in the expansion parameter M n /mN. It is 
this observation that we will use in the power counting 
as outlined in the next section. Note, if we use the ex- 
pression a = — AM^ci for the irN sigma term in the ab- 
sence of strangeness, y = 2(p\ss\p) / (p\uu+dd\p) = (with 
\p) the proton state), the expression for g\ given above 
agrees to that derived in Ref. [18] up to the replacement 
M~ — > 4M|- — c.f. appendix IA.1I The numerical values 
for gf from this derivation as well as an alternative one 
using SU(3) ChPT (c.f. appendix [Pj) agree with eq. (O 
within two standard deviations. 



3 Power counting 

It is crucial for this study to identify a power counting 
that allows a comparison of the contributions to the nu- 
clear EDMs from CP-odd transition currents to those from 
the CP-odd NN potential. The power counting originally 
proposed by Weinberg for nuclear matrix elements |51j , in 
spite of its many successful applications, is not able to ex- 
plain analogous ratios studied numerically in Ref. |52j - 
we will therefore modify it slightly, as explained below. An 
alternative scheme, which for the diagrams studied here 
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Fig. 3. Contributions to the CP-violating two-nucleon potential: (a) LO contributions, (b)-(f) NLO and N 2 LO contributions, 
where the former class contains the and the latter the g{ coupling. Solid lines denote nucleons and dashed lines denote pions. 
The CP-violating vertex is depicted by a black box. For each class of diagrams only one representative is shown. 



gives consistent results, is presented in Ref. [22]. However, 
it remains to be seen, if they are indeed equivalent. 

In Weinberg's counting, contributions to the deuteron 
EDM that come from a CP violating potential (c.f. ng.[l} 
are regarded as reducible, while the transition currents are 
counted as irreducible. Thus, one needs to power-count 
the nuclear wave functions and the photon couplings sep- 
arately, making it necessary to assign a scale to a discon- 
nected nucleon line. For dimensional reasons the corre- 
sponding 5^ function is identified with 1/p 3 , where p de- 
notes the typical momentum appearing in the evaluation 
of the integrals, identified with the pion mass, M ff . How- 
ever, if indeed nucleon momenta are of order M n , the two- 
nucleon intermediate state appearing between the photon 
coupling and the CP-violating TV TV potential is off-shell. 
Thus, also this contribution is to be regarded as irreducible 
with the two-nucleon propagator counted as (p 2 /mj\r)~ 1 , 
where denotes the nucleon mass. Again p is identified 
with M w . This power counting properly explains the nu- 
merical observations of Ref. [52] and will be used in this 
work as well. For more details we refer to Ref. [55] . 

3.1 Power counting for the contributions of the 
single-nucleon EDMs 

In a world where CP violation beyond the SM is driven 
by the 6*-term, within the effective field theory the single- 
nucleon EDMs start at the one-loop level. At the same or- 
der there are two counter terms — the hi terms in eq. (JTJ). 
The isospin structure of the loops gives that the isoscalar 
component of the single-nucleon EDMs is suppressed by 
one order in the counting compared to the isovector one [H 
[5]. However, this suppression is not present for the counter 
terms [501131] . and therefore for the power counting we 
may estimate both the contribution from the do as well as 
the di term from the estimate for the leading loop contri- 
bution given by 

g e x (M n /F n ) x (eAU) x (1/M*) x (M^f/{^) 2 
~ eg e F^M^/m 2 N , 



where the dimension-full factors in the first line come from 
the regular ttNN vertex, the photon-pion vertex (with e 
the unit of charge), the propagators and the integration 
measure, respectively, and we identified (AirFj,-) ~ mjv- In 
order to derive from this the total transition current we 
need to multiply the estimate with (1/F 2 ) x tun/M 2 from 
the TV TV potential and the two-nucleon propagator, respec- 
tively. We therefore find an estimate for the leading con- 
tribution to the total transition current from the single- 
nucleon EDM of the order of eg^j '(F n mNMn). Thus, the 
single-nucleon EDMs start to contribute to the deuteron 
EDM at NLO, as we will outline in the next subsections 
- c.f. tabled] 

3.2 Power counting of the irreducible CP-odd TV TV 
potential 

The leading diagrams for the irreducible CP-odd TV TV po- 
tential are shown in fig. [3] The leading, isospin-conserving, 
CP-odd one-pion exchange can be estimated as g\j '(M ff F ff ). 
However, as will be discussed in the next section, this term 
docs not contribute to the deuteron EDM due to selection 
rules. The first non- vanishing contribution comes from the 
subleading, isospin- and CP-odd coupling g\. It is esti- 
mated to contribute as gf/(M 7r F n ) ~ ^/(mjv-FV)) where 
we used the empirical relation, presented in the previous 
section, gf/g^ ~ M^/rtiN- This contribution will be called 
leading order (LO). 

A CP-odd pion exchange potential from a g§ coupling 
on one vertex and an isospin-odd, CP-conserving coupling 
on the other also leads to a non-vanishing contribution 
to the deuteron EDM [T51I22] . As long as we focus only 
on contributions to the deuteron EDM, the impact of 
the resulting potential is effectively a redefinition g\ — > 
g{[l + g e j3 1 /{2g A g 6 l )} [22], where /3i is the strength pa- 
rameter of the isospin-odd, CP-even ttNN vertex and <?a 
is the axial-vector coupling constant of the nucleon. The 
Nijmegen partial- wave analysis provides |/3i| < 10~ 2 [54], 
which is consistent with estimating its value from the same 
mechanism used above, namely via fj-ir° mixing — see 
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Table 1. Power-counting scales of the CP-violating AW potentials (left) and (total) transition currents (right) relevant for the 
two-nucleon contribution to the #-term- induced EDM of the deuteron. Note that the equivalence AttF^ ~ tun is assumed. 





NN potential 


(total) transition current 


LO 

NLO 

N 2 LO 


gl/{m N F*) ~ gi/(M n F n ) 
glM^/{m 2 N F^ ~ g(/(m N F v ) 
g e Mi/{m%F n ) ~ glM*/(m 2 N F*) 


g 6 e/(M%F w ) ~ g{em N /{M^F*) 
g^e/(M n m N F n ) ~ 5l e/(M 2 F T ) 
g$e/(m 2 N F n ) ~ gfe/(M n m N F n ) 



fig. |H Thus the inclusion of /?i shifts by a few percent 
at most and can therefore be neglected, given the signifi- 
cant uncertainty of g\. 



In summary, to the order we are working, the only 
contribution to the CP-odd NN potential that needs to 
be considered for the deuteron EDM is the isospin-odd 
tree-level contribution proportional to g\ and its vertex 
corrections. 



Fig. 4. Isospin-odd CP-conserving nNN vertex (black circle) 
induced by 7r°-7)-mixing and the CP-conserving r)NN coupling 
(open circle). 

The first rclativistic correction is the recoil correction 
to the qa vertex, given by — <?a/ {^itt-nF^S ■ (pi +p2)v-kr a 
where pi t 2 are the nucleon momenta and k is the outgoing 
pion momentum. It is suppressed by two orders compared 
to the gA vertex due to the additional energy dependence 
(since v = (1,0) and k = p\ — pz). 

To one-loop order there are a couple of diagrams as 
shown in fig. [3] The power counting gives for these dia- 
grams g%M^ I '(mfjFir), where we identified AttF^ ~ tun- 
Thus, the loop contributions with the CP-violation in- 
duced via the coupling g^ are suppressed relative to the 
leading, non-vanishing contribution to the potential (pro- 
portional to gf) by one power of M^/itin and therefore 
contribute to NLO. However, as outlined below, the spin- 
isospin structure of all these diagrams is such that they 
do not contribute to the deuteron EDM. At N 2 L0 the 
same topologies appear, however, with g^ replaced by g\. 
In addition, also triangle topologies of type (d) with the 

ttkNN vertex from C\jq [35] as well as vertex corrections 
(diagrams (e) and (/)) formally appear at this order. As 
shown below, besides the latter class none of the men- 
tioned diagrams contributes to the deuteron EDM. 

On dimensional grounds CP-odd four-nuclcon opera- 
tors start to contribute at order M^/m^ relative to the 
leading term. Their largest ^-term-induced contributions 
are isospin conserving (c.f. second line of eq. {!])). Thus, 
as a consequence of the Pauli-Principle, they change the 
two-nucleon spin. Therefore they do not contribute to the 
deuteron EDM. Their isospin-violating counter parts (c.f. 
third line of cq. (|T])) however, contribute, but have a rela- 
tive suppression of order (M^/tun) 2 and are therefore of 
N 3 L0. 



3.3 Power counting of the irreducible transition 
currents 

We now turn to the transition currents. As explained in 
the beginning of this section, in order to compare the 
contribution from the CP-odd NN potential to that of 
the CP-odd transition currents, the former needs to be 
multiplied by emjv/M 2 . Thus, the leading order contri- 
bution of the total transition current is estimated to scale 
as g/°em N /(M*F„) ~ g e e/(M*F*). 

The tree-level contribution, shown in fig. [SJ is formally 
of NLO, however, turns out to be of isovector character 
and thus does not add to the deuteron EDM. 

The one-loop contributions to the irreducible transi- 
tion current are estimated as c/q e/ (mj^F^) and are there- 
fore of N 2 L0. The naive power counting of the diagram 
classes depicted in fig. [5j((i) and fig. [5J(e) is slightly more 
subtle due to the cancellation of one of the nucleon propa- 
gators by the energy dependence of the tttt^- vertex. There- 
fore these diagrams are part of the irreducible transition 
current and appear at N 2 L0. 

Formally at the same order there are in addition con- 
tributions from CP-violating tree-level currents as well as 
contributions from CP-conserving transition currents scal- 
ing as (7Q<5m rap /(mArM 2 F 7r ) [52], where Sm np denotes the 
total neutron-proton mass differencial. 

In table[T]the power-counting scales of the CP-violating 
irreducible NN potentials and those of the irreducible as 
well as of the total transition currents can be found. This 
completes the discussion of the power counting. In the 
next section the various diagrams are discussed explicitly. 



4 EDMs from the term 

The computation of the two-nucleon contributions to the 
deuteron EDM is most efficiently performed in the Breit 
frame defined by q = P — P' = (0, P — P') where P and 
P' denote the total four-momenta of the incoming and 

3 We would like to thank J. de Vries, U. van Kolck and 
R. G. E. Timmermans for drawing our attention to these cur- 
rents. 
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Fig. 5. Contributions to the CP-violating transition current: (a) NLO contribution, (b)-(l) N 2 LO contributions. Solid lines 
denote nucleons and dashed line denote pions. The CP-violating vertex is depicted by a black box, a CP-conserving, but 
isospin-violating vertex by a filled circle. For each class of diagrams only one representative is shown. 



outgoing deuteron states and q the momentum of the ex- 
ternal 'Coulomb-like' photon. The electric dipole moment 
d of the deuteron is then defined by 



lim 



F 3 (q 2 ) 



2m 



N 



(13) 



where the electric dipole form factor F3 is related to the 
P- and T-violating transition current operator ( Jpj al ) M by 



(J =1,J' Z = ±1;P' |(J$ al )°| J = hJz = ±1',P) 



2m 



N 



(14) 



where J is the total angular momentum of the deuteron 
and J z and J' z its z-componcnts for the in- and out-state, 
respectively. 

The total CP-violating transition current J^ al can be 
separated into two contributions of different topology (see 
fig. [I}: two-nucleon-reducible transition currents where the 
P- and T-violation is induced by a CP-violating two- 
nucleon potential on the one hand, and irreducible CP- 
violating transition currents on the other. These will now 
be discussed in detail. 



4.1 Contributions from the CP-odd NN potential to 
the deuteron EDM 

In order for a P- and T-violating two-nucleon potential to 
contribute in the deuteron channel, it must induce 3 Si- 
3 D\ — > 3 Pi transitions, i.e. isospin to isospin 1 and spin 
1 to spin 1 transitions since the photon-nucleon coupling 
is spin independent — it therefore must be antisymmetric 
in isospin space and symmetric in spin space. 

Contributions to the CP-violating two-nucleon poten- 
tial can be further separated into irreducible and reducible 
potentials. The latter class consists of a CP-violating po- 
tential and of multiple insertions of the NN potential in 
the 3 Si- 3 Di state and/or in the intermediate 3 Pi state, 
which can be either absorbed into the deuteron wave func- 
tions, or into the intermediate NN interactions in the 3 Pi 
state and therefore do not need to be considered sepa- 
rately. 

The leading contribution to the CP-violating two-nu- 
cleon potential is the class of tree-level diagrams depicted 
in fig. |3](a). The tree-level potential induced by the g$ 
vertex is given by [T71HU] 



Ha)(0 



: 9o9A >_ 

' 2F n P + Ml 



(o-(1)-<T(2))T(i)-T( 2 ), (15) 



where I denotes the pion momentum running from nucleon 
1 to nucleon 2. It is spin antisymmetric and isospin sym- 
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Table 2. Leading order contributions to the deuteron EDM from the gf vertex without (d PW , PW: plane wave) and with (d e AIS , 
MS: multiple scattering) intermediate 3 Pi-interactions and the total leading order contribution d^_ in units of (<?i/<?o)G° efm 
with G% = gggArriN / - calculated in Zero-Range-Approximation (ZRA), with the Argonne Vis [55], Reid93 [56] and CD- 
Bonn [57] potentials. 



Potential 3 Di-adm. d e PW d e MS dl 



[T7l[2T] 


ZRA 




-1.8- 


10" 2 






-1.8 ■ 


10" 2 


nasi] 


A wis 


5.76% 










-1.43 


10~ 2 


[TSPB] 


Reid93 


5.7% 










-1.45 


10" 2 


[23P4] 


Reid93 


5.7% 


-1.93 


■ 10~ 2 


0.40 


10~ 2 


-1.53 


10~ 2 


This work 


CD-Bonn 


4.8% 


-1.95 


■ 10~ 2 


0.44 


10~ 2 


-1.52 


10~ 2 



metric and does not induce 3 5i- 3 Pi 3 Pi transitions [TjJ 

mm- 

The potential induced by the g\ vertex reads 



11 a 



L0 (l) 



■ 9i9A 



4F n I 2 + Ml 



(cr (1) +cr (2) )(r ( 3 1) - r ( 3 2) ) 



+ (t(i)-o-( 2 ))(t ( 3 1) +t ( 3 2) ) (16) 



with I as above. It is the same as in [TUHlTllTS] with g\ re- 
placed by gi . This potential-operator has a spin-symmetric 
and isospin-antisymmetric component and thus contribu- 
tes to the transition current in the deuteron channel. In or- 
der to evaluate its contribution to the EDM of the deuteron 
we resort to the parametrization of the deuteron wave 
function of [57] with a 3 P>i-state probability of 4.8%. In 
order to include the NN interactions in the intermediate 
3 Pi-state we use the separable rank-2 representation of 
the Paris nucleon-nucleon potential of Ref. [59] (PEST). 
The resulting contributions to the deuteron EDM listed 
in table H] are in agreement with the results for gi of 
Ref. [H] using the Argonne v 13 potential, of ref. [TW2"3"ll2"3] 
using the Reid93 potential and of ref. fT7Il2"l"] . where the 
deuteron wave function has been used in the Zero-Range- 
Approximation (ZRA). The 3 P>i-admixturc is found to 
enhance the deuteron EDM by about 20%, whereas the 
interaction in the intermediate 3 Pi-state reduces the con- 
tribution by about the same amount. 

Loops formally start to contribute at NLO. The re- 
ducible component of the box potential of fig. [3](6) consti- 
tutes a static one-pion exchange and is already accounted 
for either by the deuteron wave functions or by the inter- 
action in the intermediate 3 Pi -state. Its irreducible com- 
ponent may be obtained by shifting the pole of one of 
the nucleon propagators into the half plane of the pole of 
the other nucleon propagator, as outlined in [6T)ll61lR)2"] : 
i/(—v ■ pi + ie) — > —i/(v ■ pi + ie). For the sum of the ir- 
reducible part of the box potential of fig. [3](&) and the 
crossed-box potential of fig. 0(c), one finds in dimensional 
regularization in d space-time dimensions 



, 9o9 3 a 



1 



16^ 2 p 3 ^(TTo VVT+e-VL 

x T(i) -T(2)(o-(i) -er (2) ) l (17) 

with £ = I 2 /(AM 2 ). Note that the divergence has been 
absorbed by a redefinition of the four-nucleon coupling 



constant C 2 (the scale [i is introduced in dimensional reg- 
ularization) 



C 2 ° C 2 



with 



o 9q9a r 



P 3 



6L 



,d-4 



16tt 2 



-[ 7B -l-ln(47r)] 



16tt 2 
(18) 



(19) 



is the Euler-Maschcroni con- 



wherc 7b = 0.577215 ■ 
stant. 

The triangular potential of fig. |3](d) gives 



^ NLO (0 



■ 9 9A 

X T(l) ■T(2)(<T( 1 ) -<T (2) ) l (20) 



where the divergence has been absorbed by a further re- 
definition of C 2 : 



w) 9o9A 



-2L 



i 



16tt 2 
(21) 

These results reproduce those of Ref. [10]. Note that all 
gff potential-operators up to one loop as well as the four- 
nucleon-vertex operators are isospin symmetric and spin 
antisymmetric and therefore vanish in the deuteron chan- 
nel. 

At IM 2 LO there are the same topologies as just dis- 
cussed, however, with the <?q vertex replaced by its isospin- 
violating counter part g\. The triangular-potential oper- 
ator fig. \3\(d) vanishes at the considered order. The class 
of the crossed-box-potential diagrams of fig. G2(c) gives 



Ad A 



l+# 



8P3 1 16^ 2 yeo+o ~ V vr+e - 



hi 



3P 



3 1 



hi 



2 16tt 2 V \M% 

( T ( 3 1) _ r (2))(°'(l) + °"(2)) 
+ ( T fl) + T ( 3 2))( (T (1) - (T (2)) 



16tt 2 



(22) 
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Resorting again to the method presented in to 
isolate the irreducible component of the box potential- 
operator fig. \5\(b), the latter is found to be the nega- 
tive of cq. (22) and to cancel the crossed-box-potential- 
opcrator fig. E](c). Therefore, contributions to the total 
CP-violating transition current induced by the CP-vio- 
lating two-nucleon one-loop potential are absent to N 2 LO 
— not only in the deuteron channel. 

The only non- vanishing N 2 LO contributions arc thus 
the vertex corrections shown in diagrams 02(e) and (/). 
The vertex correction on the CP-conscrving vertex is read- 
ily accounted for, since we use the physical irNN cou- 
pling constant in our calculations. The situation is some- 
what different for diagram 0(e), where the physical value 
of the coupling constant is not known, but was calcu- 
lated/estimated in Sec. [2] Since <?q only appears at the 
one-loop level in the case of the deuteron EDM, we only 
need to consider g\ here. The quoted uncertainty for g\ is 
of the order of 50%. On the other hand, the correspond- 
ing correction for the CP-conserving irNN coupling con- 
stant, the so called Goldberger-Trciman discrepancy, is 
very small [63], such that we may safely assume that the 
uncertainty given for g\ is sufficiently large such that it 
includes vertex corrections. 

Thus, the only piece of the NN potential that is CP 
odd and contributes to the deuteron EDM is the tree- 
level diagram depicted in fig. [3](a), with the g\ coupling 
employed in the CP-odd irNN vertex: it is the LO poten- 
tial. 



with £ = \p' 2 — P2 + q | 2 /(4M 2 ) in terms of the initial 
(final) momentum Pi(p'i) of nucleon i and the momen- 
tum of the out-going photon q. Although the operator 
(|23[) contains an isospin-symmctric component, it is spin- 
antisymmetric and vanishes in the deuteron channel. 

The diagram classes depicted in fig. \5\(d) and (e) van- 
ish in the deuteron channel, since they are isovectors. 

In addition there are diagrams at N 2 LO where the 
photon couples to a vertex correction (fig. [5](i) and (j)); 
however, terms that contain the g$ vertex turn out to 
be isovectors and thus do not contribute to the deuteron 
channel, and those that contain g\ start to contribute only 
at N 3 LO. 

The triangular diagrams depicted in fig. 0(6), fig. 0(c) 
and fig. 03(/) arc all of order IM 2 L0. Diagrams of the types 
of fig. [5](c) and fig. 03(/) vanish in the deuteron channel 
which can be readily seen from their isospin components: 
diagram fig. 03(c) is proportional to t 3 2 ^ (photon coupling 
to nucleon 2 ) and diagram fig. 03(/) is proportional to 
2t( 2 ) + i{ T (i) x T (2))- A class of currents that has a spin- 
and isospin-symmctric component is depicted in fig.[S](6): 

Stb) ) = V ( T (l) • T (2) - T (2) 

x (Ifa-pi) (pi- Pl ) • * m + (1 2)) (24) 

with 1(1) = - arctan(|Z \/(2M n ))/(8n\l |) [35]. Resorting to 
the CD-Bonn wave function of the deuteron as used above, 
the resulting go-contribution to the deuteron EDM for the 
3 S\ state and 3 D\ admixture is found to be 



4.2 Contributions from the CP-odd irreducible NN 
transition current 

In order for an irreducible transition current not to vanish 
in the deuteron channel, it has to induce 3 Si- 3 Di — > 3 Si- 
3 Di (isospin to isospin and spin 1 to spin 1) transitions. 
It therefore needs to be an isoscalar operator, symmetric in 
spin space. Therefore the tree-level transition currents — 
c.f. fig. [5](a) — that are all isovector in character, do not 
contribute to the deuteron EDM. The relevant CP-odd 
irreducible one-loop NN current operators are listed in 
fig. 03(6)-(j). Diagrams involving CP-even JViWy-vertices 
have been neglected here since, to the order we are work- 
ing, they do not yield EDM contributions: according to 
eq. (TT4]) EDM contributions are extracted from the 0th- 
component of matrix elements of transition currents. The 
leading order, CP-even NNnj vertex ie(gA/ F%) £• Se a3b T b 
(see appendix A of [33] where 7 is the "Coulomb photon" , 
£ = (1,0)) does not have a non-vanishing Oth-component 
for S = (0,er/2). 

The diagram classes depicted in fig-03(<7) and fig-OS(^) 
are of order g® e / '(m 2 N F^) and thus N 2 LO. For a photon 
coupling to nucleon 2 the two-nucleon-irreducible compo- 
nent of diagram fig. 03(ff) an d diagram fig. f5](/i) give 



z9Wa 



l28irF 3 M n 



1 



— — arctan \/f 



(T(i) ■ T(2)-t ( 3 2) )(<t (1) -ct (2) ) . (p' 2 -p 2 + q) 



(23) 



dL. n = -2.00 • 10~ 4 xGjefm- 0.53 ■ 10~ 4 x G° e fm 
» ., ' * .. ' 



3 Si 



3 Di-adm. 



(25) 

where G° := glg A m N /F.„. 

The N 2 LO class of diagrams depicted in fig.[5J(fc) gives 

E2 



5t» ) =l 2K (T W 



T (2) ~ T (1) T (2) 



q-(i) • (P1-P1) + q-(2) • (gg - P2) 
[(px-ptf+M*] [(p 2 -p 2 ) 2 +M 2 ] 



(26) 



where Sm np is the total neutron-proton mass difference. 
It enters here because the Oth-component of the 'jtttt ver- 
tex reads Sm np (-5m np ) for an initial np (pn) state, effec- 
tively inducing an isospin- violating irir'y vertex^. However, 
as mentioned above Sm^p 6 = 2.6 MeV is already smaller 
than the naive-dimensional-analysis value eM 2 /m^ which 
is the main reason why g$ contributions are being pro- 
moted by one power in M^/tun compared to gf in the 
power counting. This is even more so for Sm np = 1.3 MeV. 
Accordingly, in order to be consistent with the rules given 
above, we need to downgrade the class of diagrams de- 
picted in fig.03(fc) by one power, to N 3 LO. The same argu- 
ment applies to the diagrams shown in fig. 03(0- The given 



4 The same effect in a different context is discussed in detail 
in ref. [44]. 
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power-counting arguments are confirmed by explicit cal- 
culations: in ref. [22] they were calculated by resorting to 
the A v\s potential and are found to be numerically small, 
i.e. 0.87- 10" 4 x G° e fm and -0.26 ■ 10~ 4 x G° e fm, respec- 
tively. In our work we take these results as an estimate for 
the higher-order contributions that are not included. 

The absence of both divergences and (undetermined) 
counter-terms up to N 2 LO ensures the predictive power of 
the two-nucleon contributions to the deuteron EDM that 
is induced by the d-term. Together with the g\ contribu- 
tion the total two-nucleon contribution to the EDM of the 
deuteron induced by the 0-term is then given by: 



= 4o 



A 6 

U N 2 LO 



-15.2 • 



9 8 o 



0.25 ±0.07 



xl0- 3 G£efm, (27) 



where the uncertainty is the one of the matrix elements 
estimated as described in the previous paragraph; the size 
of the uncertainty is consistent with what is derived from a 
naive dimensional analysis. Alternatively we may express 
the result directly in terms of 8, the strength of the QCD 
#-tcrm, and write 



'LO 



N^LO 



= - ((6.9 ± 3.7) - (0.6 ± 0.3)) x 10~ 4 9 e fm , (28) 

where the uncertainties now contain, in addition to the 
one given in eq. (|27p . also the uncertainties in the cou- 
pling constants g$ and g\. Therefore the final result is 
completely dominated by the contribution from the CP- 
and isospin-violating tree-level potential proportional to 



5 The quark-color term and the deuteron 
EDM 

The quark-color EDM (qCEDM), a CP-violating gluon- 
quark coupling, has both an isoscalar as well as an isovec- 
tor component |64j . Thus, in addition to all operators 
considered in the previous sections arising from the nu- 
cleon bilinear proportional to go in eq. ([1]), CP-odd po- 
tentials and transition currents arising from the nucleon 
bilinear proportional to g\ in eq. |T]) are now expected 
to be on the same footing. In particular, the leading- 
order contribution is now defined by the CP-odd poten- 
tial from the one-pion exchange with one c/i-vertex and 
is given in table [2] with g\ replaced by g-y. Note that 
therefore IMLO q cEDM here corresponds to LO in table [TJ 
Contributions at the one-loop level arising from CP-odd 
potential-operators are now of order N 2 LO q cEDM com- 
pared to the leading g\ potential and contributions from 
irreducible current operators are IM 3 LO q cEDM- Further- 
more, four-nucleon vertices associated with the isovecto- 
rial component of the qCEDM, which do not vanish in the 
deuteron channel, occur at N 2 LO q cEDM- It is thus suffi- 
cient to consider only those one-loop contributions with 
the largest to-be-expected yield in order to ensure that no 



extraordinary enhancements occur and the uncertainty es- 
timates are reliable. 

The relevant contributions that appear for the qCEDM 
for the CP-violating potential were already discussed in 
the previous section. The share of the deuteron EDM gen- 
erated by the class of diagrams depicted in fig. [5](6) is 
found to be equal to the one by the class shown in fig.[5](c) 
and non-vanishing. Also the contribution of the diagram 
of fig. [5](e) is non- vanishing. In addition, all these appear 
to be enhanced by a factor of 7r that commonly appears 
for triangle topologies. The other diagrams cither vanish 
(fig. O (d) and (/)) or are suppressed (fig. [5] (g), (h), (i) 
and (J)) because the ^-enhancement is missing. The nu- 
merics shows that no extraordinary enhancement of one 
particular diagram class is encountered. This confirms the 
dominance of the leading-order q\ contribution from fig. 
0(a) — see ref. [22]. 



6 Summary and conclusions 

As already stated in the introduction, the established re- 
lation between the QCD 8-teim and the CP-odd ttNN 
coupling constant is not sufficient to predict the size of 
the electric dipole moment of a single nucleon (neutron or 
proton) with the help of effective field theory, since the 
calculable one-loop contributions are of the same order 
as undetermined counter terms. However, this unpleasant 
feature is not present for the two-nucleon contributions 
of the deuteron and other light nuclei, which contribute 
already at tree-level order — unaffected by any counter 
terms — and which can be derived — admittedly with a 
large uncertainty — up-to-and-including the order N 2 LO, 
see eqs. (27]) and {28]) at the end of sect. 0T1 The N 2 LO 
contributions of these results are (up to vertex corrections 
discussed in Sec. [3]) solely governed by the irreducible tran- 
sition currents. The latter include loops which for the first 
time have been calculated in the present work. Note that 
any contribution with unknown coefficients can only show 
up at N 3 LO. 

The dominant part of the deuteron's two-nucleon EDM 
from the QCD 0-term resulted from an isospin-violating, 
CP-odd ttNN coupling constant. The isospin violation 
was estimated via fj-ir° mixing (where fj stands for the 
SU(2) or U(2) remnant of the rj and rj') and should lead to 
a nominally small contribution of the order of eM^/m 2 Nl 
if the relation between the CP-odd r/NN and tt°NN cou- 
plings had followed naive dimensional analysis. However, 
the CP-odd, isospin-conserving irNN coupling induced 
from the #-term, <7q, is found to be exceptionably small, 
such that the isospin-violating coupling gf — as already 
observed by Lebedev et al. [TH] — is effectively only sup- 
pressed by one power in the counting. 

This is important since the one-pion exchange with 
one g$ vertex cannot contribute to the two-nucleon part of 
the deuteron EDM because of isospin selection. This was 
summarized in the folklore that the deuteron would be 
blind to the two-nucleon contributions generated by the 
6*-term. This folklore, however, should be abandoned. A 
measurement of a non- vanishing neutron, a non- vanishing 
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proton and a non-vanishing deuteron EDM would suffice 
to determine the strength of the QCD 0-term, 9, from 
data. Note that the two-nucleon part of the deuteron EDM 
given in ([28)1 is in fact of the same magnitude and therefore 
comparable in size with the non-analytic isovector part of 
the nucleon EDM as calculated in ref. [31], which is, using 
as input the value of g$ from cq. ([8]), 

^non-analyt. = ^ ± g) ^ 1Q -4 q £ fm ^ (2Q) 

where the uncertainty contains both the variation of the 
loop scale as proposed in ref. [3T] as well as the uncertainty 
in <7q. This number may presumably be taken as a scale 
which governs the single nucleon EDMs. Note, however, 
that the non-analytic contribution to the isoscalar part of 
the nucleon EDM is an order of magnitude smaller due to 
a suppression by a factor M^/tun as well as the absence of 
a chiral logarithm. Whether the proton or neutron EDM 
arc really of the same magnitude as the two-nucleon part 
of the deuteron EDM is a question which only experiments 
might eventually be able to answer. 

Fact is that, under the assumption that the electric 
dipole moments are driven by the CP violation that is 
induced by the QCD (9-term, we now can give a relation 
between the total EDMs of the deuteron, the neutron and 
the proton and the calculated two-nucleon EDM part of 
the deuteron: 
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A Alternative derivations of the 0-term 
induced nNN couplings 

A.l An update of the derivation of Lebedev et al. |18| 

In addition to the usual parametrization of the 0-term- 
induced isospin-conserving and CP-violating ttNN cou- 
pling 

g e = ^(N\uu-dd\N) (A.l) 

r it 

the authors of ref. |18) introduced — via the 7r°— 77 mixing 
— the isospin-breaking counter part 



d D = d„+d p -((6.9±3.7)-(0.6±0.3))xl0~ 4 6»efm. (30) 

A cross-check of the so-extracted 9 value would be possible 

— still solely from data — by a measurement of the EDM 
of 3 He. Another strategy to test or falsify the 9 value would 
involve lattice QCD calculations and just two successful 
EDM measurements, namely one single-nuclcon EDM, i.e. 
the one of the neutron or proton, and the deuteron EDM. 
If even all three of them are measured, then one could 
use lattice QCD for a first test correlating the proton and 
neutron EDM results in terms of the parameter 9 and to 
use formula (|3"0|) for an additional, orthogonal test. 

If indeed the QCD 0-term would have failed these tests 

— cither by a direct comparison of data or by the addi- 
tional involvement of lattice QCD — then the following 
picture would emerge: in case djj — d n — d p is sizable com- 
pared to what eq. \30\) in combination with experimental 
or lattice data predicts, then the dimensional analysis re- 
veals a dominance of the quark-color EDM, feeding the 
coupling proportional to g\ . Note that ref. [22J stated the 
dominance of the quark-color mechanism already under 
the assumption that dr> — d n — d p itself is sizable. The dif- 
ference emerges since the term proportional to g\ was not 
considered in that work. On the other hand, if this differ- 
ence is very small, most probably neither the #-term nor 
the quark-color EDM is at work, but one or several of the 
other dimension six CP- violating operators [22]. More in- 
sight can be gained from a study of the EDM for 3 He. This 
reasoning stresses once more the need for high-precision 
measurements, not only of the neutron EDM but also of 
the EDMs for light ions like proton, deuteron and 3 Hc. 



A = -J- ^-rf -^(Nluu + dd-2-sslN) . (A.2) 

Inserting the strong-interaction contribution to the neutron- 
proton mass difference (m u — m,d)(N\uu — dd\N) = Sm^ 6 
and utilizing the parameter e as defined at the beginning 
of sect.[2J we derive (J5J again: 



9o 



6m%*(l - e 2 ) 
4-F 7i-£ 



Similarly, starting now from eq. (|A.2[) . we get 



a - ~ 6 (l A) , M * 

9 '~8F n {1 e > e M*-M* 



m(N\uu+dd-2ss\N) (A.3) 



with Ml = 2Bm+0{M 2 ) and M 2 K = B{m s +rh)+0{M 2 ) 
for the square of the pion and kaon mass, respectively, 
where here M. is the quark mass matrix for three light 
flavors. According to refs. [^TllrJrj] we have rh(N\uu+dd— 
2ss\N) = rh{N\uu + dd\N){l - y) with ovzv = (T„ N {Q) = 
m(p\uu + dd\p) and y = 2{p\ss\p) / (p\uu + dd\p) , where \p) 
denotes here the proton state. The final result is therefore 



91 ~ 8^ (1 €)e M 2 K -M 2 



a wN {l-y). (A.4) 



5 Actually, via the 7r -?7g mixing. For consistency, we re- 
placed here their 7r°-?7 mixing angle by the customary one of 
chiral perturbation theory |65] — note the explicit rh subtrac- 
tion in the denominator. 
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Table 3. The value of go, gf, and the ratio gf/go predicted from eqs. (|A,9|) and (|A.10|) with (i) the original SU(3) parameters 
bo and &f of ref. [35], with (ii) the alternative set of parameters based on eqs. (|A.13[I and (|A.14|) . (iii) in the case that bo + b F 
of (i) are replaced by C5 of eq. (|A.12[) . The listed uncertainties do not contain systematical SU(3) errors. 









g e [eGeV- 1 ] 


gt [^GeV- 1 ] 


B 1 9 

9i/ go 


(i) 


br> & b F 


from |35| 


-0.026 ± 0.002 


0.00092 ± 0.00017 


-0.036 ±0.007 


(ii) 


bo & b F 


alternative 


-0.023 ±0.005 


0.00088 ± 0.00016 


-0.038 ±0.011 


(«i) 


b D +b F - 


-> cs 


-0.018 ±0.007 


0.00092 ± 0.00017 


-0.051 ±0.022 



Inserting 5m s ^ s = (2.6 ± 0.5) MeV from ref. [42], F v 
92.2 McV, and the ITS quark masses at 2 GeV from [3< 
we get 



and 



gl w (-0.018 ± 0.007)6* 



g{ « (0.0012 ±0.0004)0 



(A.5) 



(A.6) 



the SU(3) counter parts of eqs. (JHJ) and (|A.4|) are 

-6f), (A.9) 



with fr ff jv(0) = 45 MeV and y = 0.21 ± 0.20 from [55] as 
additional input. 
Thus we find 



ii. 
9% 



Ml 



2 Ml-Ml 



^j v(0)(1 ~ y) 

jC Strg 



-0.07 ±0.04 



(A.7) 

as the ratio of the isospin-brcaking versus the isospin- 
conserving CP-violating ttNN coupling constants which 
are induced by the 6*-term. If we rather applied the val- 
ues (7^(0) = 59(7) MeV and y w from ref. [57] (for an 
update of this work see ref. [55]), we would get 



g{ « (0.0021 ± 0.0004) and 



9 6 o 



-0.11 ±0.05 
(A.8) 

as values for gf and the ratio instead. Note that the ratios 
listed in (|A.7|) and (|A.8|) arc compatible with the estimate 

m 



9n°NN 



AB6m* 

Fir 



and 



9n°NN 



ABOm* 
~FZ 



respectively, whereas 
-AB6m* b d 



9^NN 



and 



9r/NN 



ABOm, 



arc the corresponding r/NN (actually ij^NN) counter parts. 
Note that AB bp and AB bp are the coefficients of the an- 
ticommutator (D-type) and commutator (F-type) term of 
the quark mass matrix with the baryon matrix. Therefore, 



gg = 4^ J+M = ^ (1 _ £2)(bg 



<h 



Fr 

AB8m* 36 F - 



b D V3 



m d -m u 



V3 



-m 



M 2 

?-^(l-e 2 ) (3b F -h 



eMl 



4(M| - Ml) 



(A.10) 



where (v / 3/4)(m £ ; — m u )/(m s — rh) is the n°-ri (actually 
TT°-r]s) mixing angle. Thus, in this case we get the ratio 



ii 

9 6 



eMl 



4(M 2 



K 



36, 



Ml) b 



>D 



D 



(A.11) 



If the values b F = -0.209 GeV" 1 and b D = 0.066 GeV -1 
of ref. [35] are inserted, we get the first row of table [3] 
Note, however, that there is a mismatch by a factor 1.5 
approximately between the SU(3) octet quantity 



b D +b F = 



ms — ms 



(-0.143 ±0.004) GeV -1 



4(M| - Ml) 

used in [5"l l3Tjll3l"] and the SU(2) low-energy coefficient 



C5 



AMle 



-0.097 ±0.034) GeV" 



(A.12) 



although according to SU(3) ChPT both quantities should 
agree to leading order, see eq. (27) of ref. 17 



A. 2 Derivation via SU(3) chiral perturbation theory 

In SU(3) ChPT the D-type and F-typc CP-violating tt NN 
coupling constants are (see e.g. the U(3) ChPT calculation 
of ref. [32]) 



Moreover, an alternative procedure to parametrize the 
above sum is 



b D + b F 



SM BtIg 



A(M 2 K+ -(M 2 + -M 2 ) 
(-0.126 ± 0.024) GeV -1 



M 2 K0 ) 



(A.13) 



where the electromagnetic mass shifts are removed (via 
the Dashen theorem [71] in the denominator) and where 
the prediction falls in-between the original one and the C5 
value. Using an analogous parametrization for bp, we get 



8{M 2 



M 2 K0 ) 



-0.196 GeV -1 
(A.14) 



6 In fact, the latter equation which is based on eq. (5.7) of 
ref. [70] predicts that the NLO correction to C5 is much larger 
than c 5 (or b D + b F ) itself, namely Ac 5 = 0.49 GeV -1 . This 
quantity is of similar size as zici = ±0.2 GeV -1 . 
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and b D = +(0.069 ± 0.024) GeV" 1 from (TOp instead of 
the above listed values from [35], such that the values in 
the second row of tableware generated instead. Note that 
the result for 2>bp — bu is approximately the same in both 
paramctrizations, namely —0.69 GeV -1 in the original one 
[55] and —0.66 GeV -1 in the modified one. 

Finally, replacing bp + &f of [35] by C5 of eq. (|A.12|> . 
we get the values in the third row of table [31 

Note that only the last SU(3) value of the ratio <?i/<?o 
is in the range of our estimate (fT2"]) , but all three are com- 
patible with the estimate of (|A.7[) . The quoted numbers 
of table [3j however, do not contain a systematical error 
connected with an SU(3) ChPT calculation. For standard 
quantities such an uncertainty is certainly of the order of 
50 %. For the quantity C5 this uncertainty should be rather 
100%— 200%, see e.g. footnote© Taking these SU(3) er- 
rors in account, the estimates of table [3] are compatible 
with the range quoted in ([12")) . 
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